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We construct the perturbative solution of the charged black hole in the Kaluza-Klein spacetime 
with the matched asymptotic expansion method. The corrections to the thermodynamic variables 
are calculated up to the post-Newtonian order. We confirmed that the method can work very well 
in the Einstein-Maxwell theory. 



I. INTRODUCTION 

In the last decade, the higher dimensional black holes have been known to have much richer structure than in four 
dimensional spacetime The discovery of the black ring solution 2, 3] revealed that the uniqueness and the horizon 
topology theorem in four dimension arc no longer valid in higher dimensions. 

The similar breakdown happens in the Kaluza-Klein spacetimes. When the mass scale is smaller than the scale of 
the extra dimension, there are the black string solutions and the black holes localized in the Kaluza-Klein spacetime 
(the caged black hole solutions) for the same total mass, which have the different horizon topologies. Furthermore, 
there exist the uniform black strings (UBS) and the non- uniform black strings (NUBS). The NUBS phase bifurcates 
from the UBS phase at the marginally stable mode of the Gregory-Laflamme instability [1, @. This system admits 
the interesting phenomena such as the topology changing phase transition from the NUBS phase to the caged black 
hole phase, and the critical dimension for the stability of the NUBS phase @. 

Since it is difficult to construct the exact solutions of the black objects with the nontrivial horizon topology, the 
numerical or perturbative approach is helpful in the analysis of higher dimensions. In the Kaluza-Klein spacetime, the 
small caged black holes have been first studied perturbatively using a single coordinate patch 0, [|| . Gorbonos and Kol 
used two coordinate patches and introduced the systematic procedure to obtain the solutions perturbatively [9[ . This 
perturbative method is called the matched asymptotic expansion (MAE). The MAE can be applied to the spacetime 
having at least two separate scales. Using the MAE, they constructed the small caged black holes up to the post- 
Newtonian order in the asymptotic zone [3, [13] ■ This method is also useful in the ultra-spinning limit of the black 
objects [Hm. 

The effective field theory (EFT) method [l3| is another useful method to systematically obtain the thermodynamic 
properties of such spacetime, which is considered to be equivalent with the MAE method. In the EFT method, the 
small scales are integrated out to give the effective Lagrangian for the coarse-grained black objects. This is applied 
to the caged black holes [TJ, [l5| . The caged black holes with charge is also studied by the EFT [r| . 

In this paper, we will generalize the MAE method to the Einstein-Maxwell theory and construct the perturbative 
solution of the charged black holes in the Kaluza-Klein spacetime. In the neutral limit, the solutions here reproduce 
the neutral results [1(| ■ We also calculated the correction to the thermodynamic variables up to the post-Newtonian 
order which confirms the EFT result 

The organization of this paper is as follows. In Sec. II, we explain the setup of the caged black hole and the 
matching method. In Sec. Ill, the Newtonian order is solved at the asymptotic zone. In Sec. IV, we solve the 
near zone perturbation. In Sec. V, we match the near zone results with the results in Sec. Ill and evaluate some 
of the thermodynamic variables. In Sec. VI, we solve the post-Newtonian equation at the asymptotic zone. In 
Sec. VII, we compute the PN correction to the global charges, and analyze the thermodynamic properties of solutions. 
We summarize our work in Sec. VIII. In Appendices, we will present some useful formulae and the details of the 
computations. 



II. MATCHED ASYMPTOTIC EXPANSION 



In this section, we introduce the method of the matched asymptotic expansion in the (n + 3)-dimcnsional Kaluza- 
Klein spacetime with the Einstein-Maxwell theory. The action is given by 

S = f(j^R-\F^)^gt«*x, (1) 
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where F^y is the field strength of the Maxwell field defined by 
unit of G = 1, hereafter. The Einstein equation becomes 
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d, L A u — d u An. In this paper we choose the 
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where the energy-momentum tensor is given by 



T, 



i 1 ' ' 



F^ X F V \ - -FapF^g^, 



and the field equation for the Maxwell field is 



0. 



(3) 
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We consider a static localized black hole solution in the Kaluza-Klein spacetime which is called "caged black hole" . 
Then we take the small black hole limit in which the scale of the black hole po is much smaller than the compactification 
scale L. Since the Maxwell charge is bounded above by the extreme limit, it can be assumed to be the same order as 
the mass ~ pQ. Therefore, we can perform the similar treatment with the neutral case [9l Hoj|. 

In this small black hole limit, we have two asymptotic zones. The spacetime asymptotes to the flat Kaluza-Klein 
spacetime in the asymptotic zone {p po), and the charged spherical black hole in the near zone (p < L). Let us 
introduce two coordinate systems. In the asymptotic zone, it is better to use the cylindrical coordinates (r, z) as 



ds z 



-dt 2 + dz 2 



dr 2 



r 2 dnl 



(5) 



in which the z-direction is periodic as z ~ z + L. The periodicity requires that we must consider the gravity from 
the mirror images of the black hole. This is equivalent to the infinite array of black holes in the higher dimensional 
Minkowski spacetime. In the near zone, on the other hand, it is better to use the spherical coordinates (p, x) m order 
to see the black hole perturbation as 



ds 2 = ~dt 2 + dp 2 + p 2 (d X 2 + sin 2 xd^l)- 
In the asymptotic zone, the relation between two coordinate systems is 

r = psinx, z = pcosx- 



(6) 



(7) 



We expand the metric and the Maxwell field in both zones. In the near zone, we consider the black hole perturbation 
with the expansion parameter 1/L as 



where h\ T 
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0(1/ L k ). gjj^ is the metric of the (n + 3)-dimensional background charged black hole spacetime 

1 



~f(p)dt 2 



Hp) 



dp 2 + ffdnz +1 , 



f(p) 
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(9) 
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and cr is a dimensionless charge parameter which becomes \a\ 
hole is 



4 p 2 ™' 

1 in the extreme limit. The total mass of this black 



Mn = 



16tt 
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where lo^ — 2n 2 /r(£±i) is the area of the k- dimensional unit sphere, S . The background Maxwell field is 



(12) 
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where the total charge of the background black hole is given by 



(13) 



and non-vanishing component of the field strength is the electric field 



/n(n+l) pS 

Ftp ~ \l^2^°^ = £ - (14) 
In the asymptotic zone, we consider the post-Newtonian expansion with the expansion parameter po as 

oo 

5^ sym) = ^ + E^ sym ' m) ' ( 15 ) 

m— 1 

where ft,( as y m > m ) ~ 0(p™) and 7/ M „ is the metric of the in + 3)-dimensional Kaluza-Klein spacetime ([5]). The Maxwell 
field is also expanded in the same way. It is worth noting that the Newtonian approximation is the order of pp because 
the mass M is written as p$ ~ GM . 

Now, it is ready to consider the matching procedure between the two zones. In the limit po "C L, there exists the 
overlap region (po *C p "C L ), in which we will match the two expansions order by order. The "matching ladder" can 
be understood by the dimensional counting. If we take the harmonic gauge, the linear part of the post-Newtonian 
equation in the order of p™ becomes 

where □ = -q^d^d^. Since we consider the static spacetime, it reduces to the Laplace equation. The homogeneous 
solutions with a multipole I behave as p l and p~ z ~™. For the asymptotic zone, the near zone black hole seems to 
have multipole moments and behave as p~ l ~ n . From the dimensional counting, the terms with p~ l ~ n in asymptotic 
solution of 0(p™) appears in the dimensionless form 

m 

h (asym, m ) Hi (17) 

jjm—l—npl+n v 7 

Therefore, the near solution of 0{1/L k ) contributes to the asymptotic solution with the multipole I of 0(p^ +n+l ). 

Similarly, the gravity from the mirror images of the black hole affects the near zone geometry in the form of multipole 
moments at the infinity, which are proportional to p l . The dimensional counting shows 

Mr r ' fe) 3 ^jt-- (is) 

Then, the asymptotic solution of 0{p™) determines to the near solution with the multipole I of 0(l/L m+l ). 



III. NEWTONIAN POTENTIAL IN ASYMPTOTIC ZONE 



We start the matching of the near zone solutions with the Newtonian order solution in the asymptotic zone. We 
will omit the script "near" or "asym" for the brevity. As explained in the previous section, the near zone background 
metric of 0(1/L ) gives the correction to the asymptotic metric in O(pq). At this order, the solution is just the 
Newtonian potential from the array of the point source as in the neutral case @. In the harmonic gauge dM 111 = 
where = — ^hr]^, the linearized Einstein equation is 

- \nh$ = StttW. (19) 

Since is written in terms of the square of and the leading order of F^ v is the Newtonian order, the Einstein 

equation is the same as the vacuum case at this order (namely, T^u = 0). The boundary condition in the overlap 
region is determined by the behavior of the near background metric. In the overlap region, the near zone background 
metric in the harmonic coordinates is given as 0(1/ L°) terms in Eq. (|D8|) . which becomes 

- -1 + 4 (20a) 

pn 
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Requiring the periodicity in z-direction and the harmonic condition, the homogeneous solution with a pole at the 
origin is given by 



oo 1 

hif = $ = p" V — — -, (21a) 

" H0 ^ (r 2 + (z-fci) 2 )t ' V ; 



A;— — c 

"ij 



4° = Z*V (21b) 



Since $ ~ Po/p n as P "C L, the coefficients are determined by the match with Eq. (|20p . This is equivalent to the 
Newtonian potential from the infinite array of point masses, 

oo 

T t [ n) = M 6 n+1 (x)6(z-kL). (22) 

k— — oo 

Imposing the Lorenz gauge d^A^' = 0, the Maxwell equation is 

= 0- (23) 

In the same way, the solution is given by 

An) = 6 = Qo 1 (OA) 

* V n^+i ^ (r 2 + (z-fci) 2 )^ V ' 

This is also equivalent to the potential obtained from the infinite array of point charges, 

oo 

Jt n) = QoS n+1 (*)5{z-kL). (25) 

k— — oo 

Now we can see that the gauge potential is proportional to the Newtonian potential, 



nUJn+lPo ' 327TO 

<E> is expanded in the overlap region for r, z <C L, 

„n n 2k 



Qo n+1 . . 

$=i/— — (26) 



* = + E 2 C(- + 2fc)|^C 2 l(cos X ), (27) 

P k=0 



where CY(x) is the Gegenbaucr polynomial (see Appendix lAl for the definition). 



IV. NEAR ZONE PERTURBATION 



We now proceed to the leading order in the near zone. The leading order correction in the near zone comes from the 
asymptotic solution of 0(pq ) in the previous section. In this paper, we only consider the monopolc (/ = 0) correction 
which is 0{l/L n ). 

Here, we calculate the linear perturbation around the black hole metric of Eq. ([9]) as g^ v = gj^ H ^ + h^ u . Since we 
consider the static perturbation, we can set the following ansatz for the metric perturbation 



ds 2 = -f{l + A )dt 2 + r 1 (l + Bo)dp 2 +p 2 {l + E )dfi 2 n+1 , 
and for the Maxwell field perturbation, 

Qo 



nu n+1 p" 



+ a ) dt 



We write the linearized Einstein equation as 



5R[iv = 8n 



ST^ - 



n+ 1 



n + 1 



(28) 



(29) 



(30) 



where T^J is the energy-momentum tensor of the background. The components of are given by 
_ t(0)« _ _ T {0)r _ T (0) X _ T (o)6 _ I 2 _ n(n + l)a 2 p 2n 



where {#i}i=i..n are coordinates of S n . The explicit forms for each components of SR^ and 8S^ V are given in Eqs. (|C9 
and (jClOp . The Maxwell equation becomes 



d p [ P n+1 a' - ^-((n + l)^o - 4, - Bo)] = 0. (32) 
zw„+i 



This can be easily integrated as 

£ 



a' = -((n + l)E -A -B + 2a ), (33) 



where ag is an integral constant. Using the above and Eq. (|C10[) , SS^ becomes 



n 2 n 2 n 2n 

6S*t = 4 Z +2 [(n + l)E -A + 2ao], (34a) 
nrr 2 n 2n 

SS% = — -4^l(n+l)Eo-Ba + 2ao] (34b) 



and 



2 In 

5S\ = -^-^[E -A -B + 2a ], (34c) 

where SS^u = g^ Bii ^ La SS a ^. Subtracting the pp-component from the tt-component in Eq. ([50)1 . we obtain the following 
equation 

pE' ' + 2E' -A' o -B o = 0, (35) 
which is the same with the neutral cases. Another independent equation comes from the xX~ com P onGn t hi Eq. (1301) . 

\/{lK + f'K) + + l)pfK + \pf{A' - B' Q ) 

+n(E - B ) = — '\[{n + 1)E - B + 2a }. (36) 
4 p zn 

Integrating Eq. (|3"5j) . we have 

pE' + E -A -B = -2C 1 , (37) 
where C\ is an integral constant. Imposing 

A a + B + (n-l)E =0 (38) 
as the residual gauge condition, Eq. (f37l> can be solved as 

Eo = ^ + (39) 
n p n 

where Ci is an integral constant. 

To solve Eq. (|36| in the gauge condition (f38j) . we introduce a new variable defined by 

* = ^(Ao + nB ) = E - B , (40) 

n — 1 

where A and B Q are written by 'J as 

A = - n£ , B =E Q - (41) 
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Using the above and Eq. (I39I) . Eq. (|36j) becomes 



p (f^Y + nf^ = 0^[a 2 (a o + C 1 )+nC 2 ] (42) 



and then the solution is given by 



2/i 



f* = ^-^l**(ao + C 1 ) + nC 2 ], (43) 



where C3 is an integral constant. Even after imposing the gauge condition (|38[) . there is still remaining gauge degree 
of freedom p — >• p + SCp^ /2p n ~ 1 which changes Eqs. (|31?|) and (|4"3")l as 

E ^^ + ^, (44a) 
n p n 

„2n 



f< n n 2 



where C 2 = C 2 + 5C and C3 = C3 + nSC. Therefore C 2 can be the pure gauge and we define the gauge invariant 
combination C3 = C3 — nC 2 . Here, C3 is the ambiguity of the bare mass and ao is the ambiguity of the bare charge 
which we can freely choose. 

As the boundary condition at the horizon, we impose the regularity conditions for Ao, Bq and Eq which fix the 
horizon position. We write the position of the outer event horizon as p+, which is the larger root of f(p) = 0, 



L»(1 + V /1- CT 2). (45) 



P + 2 

The regularity of at p = p + is guranteed, this condition is satisfied if the right-hand side of Eq. (|43| vanishes as 
p — > p+, which gives 

-<j 2 (a + C 1 ) + -nC 2 = (C 3 +nC 2 )?±. (46) 
ll p 

In general, the perturbed spacetime may not become extremal or neutral when the background spacetime is extremal 
[a = 1) or neutral (cr = 0). Then, we will further fix the parameter region of a so that the solution becomes extremal 
when cr = l and neutral when a — 0. The extremal condition is attained by imposing that the right-hand side of 
Eq. (|43p has the double roots, (/^ , )'|p =p+ . cr= i = 0, which becomes 

(a + Cx + nC 2 )\ a=1 = = (C 3 + nC 2 )\ a=u (47) 

where we used Eq. (|46[) . The neutral condition is 

(2C 3 + nC 2 )\ a=0 -0. (48) 

Although Ci and ao may depend on a in general, we can set these to satisfy the above conditions for not only a = 0, 1 
but also arbitrary a as 

C 2 = C 3 = 0, a = -C x . (49) 

Therefore, "J becomes to be zero, that is, 

* = 0. (50) 
Then, from Eqs. ([3U1) . (pJ9"| and (|4"Tj) . the perturbation becomes 

(51a) 
(51b) 

(51c) 

- + ax, (51d) 

where ai is an integral constant. As a result, the near solutions depend on two parameters Cx and ax, which are 
determined by matching with the asymptotic solutions in the overlap region. 



Ao = 


-2Cx 


B = 


2Cx 




n 


E = 


2Cx 




n 




a 


a Q = 
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V. MATCHING FROM THE NEWTONIAN POTENTIAL 

As mentioned in Sec. UH the near zone monopole correction of 0(1/L n ) comes from the monopole parts of the 
Newtonian order asymptotic solution. The monopole from the asymptotic zone behaves as a constant which has the 
dimensionlcss form, po/L n at the Newtonian order. 

Since the leading order terms in the near zone perturbation are not affected by the gauge transformation into the 
harmonic coordinates (|D7|) . the near solution in the overlap region (p <C p) becomes 

h\r T ' n) = ~fA ^2C u (52a) 
hfr' n) = ^ - — , (52b) 

^(near.n) = 2^ = ^ 2 (52c) 

n 

On the other hand, the p^/L n term in the expansion of the Newtonian potential (|27p is 

hlr ym ' n) = t> = ^ + 2«n)<^ + ... (53) 

and then the matching provides 

Ci=CWf|- (54) 

We now define the expansion parameter A = C( n )Po /L n . Also, for the Maxwell field the leading order term of the 
near solution matches the gauge potential <p of the asymptotic solution in the overlap region as 

a (ncar,n 1 

A ~ ai = -A. (55) 

nLu n+1 p% 

Then, combining Eq. (|5"Tj) with the ansatz (|28|) and ([29]) , the near zone metric up to 0(1/L n ) becomes 

sir = -/(1-2A) + 0(1/L 2 ™), (56a) 

9 { ;r i} = + ^(Vi 2 "), (56b) 



g^ = p 2 ( 1 + — ) + 0(1/L 2n ), 



and the gauge potential in the near zone is 

A (nca, = + JXQ^ 0(1/L 2 "). (57) 

nuj n+ ip n nui n+ ip% \ 2p n J 

From the leading order near solution (|56[) . the correction to the local constants can be computed. The surface gravity 
k becomes 

1 \d P 9tt\ 



2 V~9ttg PP 
2p\ +lW \ 

= k ( 1 - — a) , (58) 



where kq = \/T— er 2 n/3g/2/3™ + is the surface gravity for L = oo. The horizon area A becomes 



A = uj n+1 (g xx ) ^ \ p=p 



= oj n+lP l +1 



-i-i 

"+ 

, n+1. 
1 + A 

n 



= Ao[l+ — A ) , (59) 
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where Aq is the horizon area for L — oo. 

We can also calculate the electrostatic potential U, 

U ^ A t \ p=p+ - A t \ p=QO = ^Ln +V /IT^ A) 

nu) n+ ip + 

= U (l+ a/1-<t 2 A), (60) 

where Uq is the electrostatic potential for L = oo. We note that the boundary condition of the asymptotic solution 
leads to A t \ p=00 = 0. 



VI. MONOPOLE MATCHING IN POST-NEWTONIAN ORDER 



The leading correction to the mass and tension in the asymptotic zone are computed through the monopole pertur- 
bation of the post-Newtonian order, 0(p^ n ). The near zone perturbation of Q(l/L n ) considered in the previous section 
gives the boundary condition in the overlap region, which behaves as p§™ /L n p n . First, we present the post-Newtonian 
equation for the gravity and the Maxwell field and then solve them. 



A. Post-Newtonian order perturbation equation 



In the asymptotic zone, we consider perturbations of the metric and Maxwell fields to the post-Newtonian order, 

g»u = V„u + ft W + , Ap = 4") + 4 2 ") , (61) 

where h)}^ and Ap are the post-Newtonian corrections. In the harmonic gauge, the post-Newtonian equation 
becomes 



(62) 



where we used the Einstein equation of the Newtonian order, G p „ = T p „ = 0. R^l is the 2nd order perturbation of 



the Ricci tensor, which gives a source term from the Newtonian order of Eqs. (|21[) and ([24]) . The Maxwell equation 
in this order is given by 



p{n)v _ u(n)va p(n) _ . (n) p(n)va 
1 H " r aft "/ia r 



0. 



(63) 



Using Eq. ((Ml), then, Eqs. (JHU) and (JM) become 

□ 



0. 



(64a) 



□ 



(2n) 



1 

2^2 



no 



2n 



(1 - ° 2 ) 



(64b) 



□ 



.4 



(2n) 



2nw„+i pJJ 



0. 



(64c) 



The general solution is constructed by the inhomogeneous and homogeneous solutions. The coefficients of the homo- 
geneous solution are determined by the boundary condition. We write the general solution as follows 



(2n) 



(2n) 



. 1 - — I $ 2 fc + Sij — $ + Pf 
2n 2 V 2 J J L" 



G(a:,a:')<SyV)d" + V , 



,(2n) _ 



2nw„+iPo 



* 2 + -A^*. 



(65a) 
(65b) 
(65c) 
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where Pf means the finite part of the integration and St, Sij, sa are the dimensionlcss coefficients of the homogeneous 
solution. The Green function in the compact space, [-L/2, L/2] x R™ +1 , is given by 

1 00 1 
G(X > X>) = ^ ((x-x0 2 + (.-.'-mLnf (66) 

where x is the coordinate vector of R n+1 , which gives r = |x|. We write the integrand as 

Si, = [(1 - * 2 )($,i®,i) - ■ (67) 



B. Matching from the near solution 



Now, we determine the coefficients, s t , s^, sa of the homogeneous term by matching with the near solution. The 
monopole moment of the near solution behaves as ~ p n - Then the relevant terms of this order in the overlap region 
have the dependence of p^ n /L n p n in the dimensionlcss form. Transforming Eqs. ([56)) and ([57)) into the harmonic 
coordinates (|D7|) . the near zone solution up to the relevant order becomes 



(near) 

9tt 



1-2A- (1 -2X)^- 
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(near) 



A 



2\ 
1 + — 

n 

2XQ 



2A\ pi 
1 A P0 



nui n+ ip nuj n+ xp 
We extract the monopole moment of 0(1/ L n ), 



+ 0(l/L 2n ,l/p 2n ), 

S tJ +0(l/L 2n ,l/p 2n ), 

-(l-X) + 0(l/L 2n ,l/p 2n ). 



n J np 

h 



(near) 

9tt 



(near) 

9ij 



2n 

2C(n) Pl n , 
n 2 L n p n y ' 

C(n)Qo P 2 n 



noj n+1 p% L n p n 

Meanwhile, the corresponding terms in the post-Newtonian solution (|65|) are 



u&n) 



(2n) 



A2n) 



S t - 2C(n) ( 1 + — 



2n 

Po 

L n p n ' 



2C(n) 



1 - 



2n 

Po 

L n p n ' 



2C(")Qol Po" 



TlUn+lPo 



L n p n 



(68a) 
(68b) 
(68c) 

(69a) 
(69b) 
(69c) 

(70a) 
(70b) 
(70c) 



The finite part of the integration term in Eq. (|65b[) does not contribute to the monopole moment in the overlap region. 
The matching between two solutions shows 



st 



$A 



a 2 C(n), 



— ((n)5ij = —S^, 
n n 



C(n)Q, 







nuj n+1 p'^ 



(71a) 
(71b) 

(71c) 
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Therefore, the post-Newtonian correction is determined as 

hT ] = -i(l + y)$ 2 + ^A$, (72a) 

A? n) = + Q° A $, (72c) 

VII. GLOBAL CHARGES AND THERMODYNAMICS 

In this section, we compute the global charges and confirm that the first law of the thermodynamics holds for the 
current cases. 



A. Global charges 

Now, we calculate the post-Newtonian correction to the global charges. The (n + 2)-dimensional mass M and 
tension r of the asymptotically Kaluza-Klein spacetime are determined by the asymptotic behavior (l7l HH , 

h [r ym) - (73) 

AST - (74) 

Then we find 

M=^(nc t -c z ), (75) 

r = ^( c *- nc -)- ( 76 ) 



The total electric charge Q is determined in the same way, 



1 



A ( asym) * -^r (77) 

(71 — 1)cj„L r n 1 



To extract the global charges, we take the limit, r ^> L,z. As in Eq. (|B1|) . $ becomes 



* , (78) 
in — 1 w n Lr" 1 



Since the Green function in Eq. has the similar behavior as <£>, the integration becomes 



G(x,xOS y (z'K + V~-- - — I / ^(z'K+V (79) 

|,|]xr»+i (n-l)w n Ir n i[_i^ ]xR n + i 



Then, Eq. ([72]) gives 



h ( asy m,2n) _ nuynQ- 2 C(n) p 2 ," 1 



(n - l)w n L n+1 r"- 1 277,(71, - l)w„ir™- 1 



-L/2 



(81) 



The second term in iS zz (a;) does not contribute to Eq. (JBTj) because it is the total derivative. The finite part of the 
integration is computed in Ref. [7i| 

Pf ( jd^xf jzA = _("-iKW^ + ipg". (82) 
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As a result, the ADM mass and tension become 



m = (" + 1 K + ^V 1+ i (1 + g2)A 

lD7r y 2 

= Mo (l + l(l + a 2 )\\ , (83) 

TL = 32^ {l ~ a )A 

= |Af (l-a 2 )A. (84) 

From Eq. (|72c[) . the total charge Q becomes 

Q = Q (1 + X). (85) 

At a = 0, these results reproduce the results in the neutral case [Io| . 

Moreover, in the extremal case [a = 1) we can confirm that our results correspond to the exact solutions in Ref. [l9j |. 
The exact solutions of the (n + 3)-dimcnsional extremal caged black hole are given by 



where 



ds 2 = -H- 2 dt 2 + H 2/,i fcdxW, A.dx" = -\ ——H^dt, (86) 

' 87T7T, 



f<*') = i+e £ (87) 

If we set p = Pq(1 + A)/2, at a = 1 the perturbative solutions constructed agree with the above solutions up to 0(A). 

There is the apparent difference between our results and the results by the EFT calculation [la ]. This simply 
comes from the ambiguity of the parameterization. For example, see Eqs. (3.7) and (3.25) in Ref. [l6[ which are 
M and r. Changing the gauge condition in Eq. (|4"9")l gives the different length scale po —> po{l + ei(c)) and charge 
parameterization a — s- er(l + £2(Aer)). We note that, to keep a = 1 to be the extreme limit, one requires £2(A, 1) = 
further. To reproduce the results by the EFT, we should impose the following gauge condition, instead of Eq. (|49j) . 

C 2 = — ?±^,C 3 = C 1 ,ao=0. (88) 

B. First law and the Smarr formula 

It is ready to confirm that our solutions satisfy the Smarr formula 

n + 1 

nM = — —kA + tL + uQU, (89) 

8-7T 



and the first law 



dM = rdL + —ndA + UdQ. (90) 
8tt 



Collecting the thermodynamic variables from Eqs. ([55)1 . (|5"9"|) . (|60l) . (|55|) and (JH3|, we write 



M = M a \ [ l + -(l + <j 2 )\), A=A \l + —^-\), Q = Q (1 + A), 
k = kq [1 — ————A J , [/ = U (l + y/l - ct 2 A), tL= ^M a (l-a 2 )X, 



(91) 



where 



(n + l)w B+ ip8 „ „+i f 1 + VI - ff 2 \ " ^ . , n(n + l), 



Mo = > 16 ; ■ -A) - ^+iPo +i (^^^ J , Qo = V 32 ^ - "n+i/>gg 

*> =^rVT^, ^ = ^V. (92) 
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From the direct computation, the above expressions are easily confirmed to follow Eq. (|89|) up to 0(A). 

To confirm the first law for the current case, we take the variation of M with L, A and Q. Changing the variables 
from (po, A, a) to (L, A., Q), the variation becomes 



{d L M) AQ 
(d A M) L , Q 
{d Q M) L , A 



d(po,X,a) 
d{L,A,Q) 



(d po M) X ,a 

(dxM) P0 . a 
(d a M) P0:X 



d(po,X,a) 
d{L,A,Q) 



nM/po 
M (1 + ct 2 )/2 
MqctX 



(93) 



Using L = C(n)»p X~ 
d(p ,X,cr) 



: - , the Jacobian d(po, A, a)/d(L, A, Q) is computed as 
d(L,A,Q) 



d(L, A, Q) 



d(p Q ,X,<r) 

PoX/L —nX/L 
Pok (1 - (2n + l)X/n)/8xnM Xk /8ttM 

poU (l - 2X)/nM XU /M U {1 + VT^a 1 )^/!^ /aM 





-ctko/8ttMo 



Using this, we can show that the variation of M satisfies the first law (|5U|) up to 0(A) 



(94) 



(95) 



Note that this is also achieved through the Harrison transformation from the neutral seed solutions [20|. On the 
other hand, our argument gives us a direct confirmation using the direct construction of the perturbative solutions. 



{d L M) A . Q 




T 


(d A M) L , Q 




k/8tt 


(d Q M) L . A 




U 



VIII. SUMMARY AND DISCUSSION 



In this paper, we constructed the perturbative solution of the small black holes with the Maxwell charge in the 
caged spacetime using the matched asymptotic expansion. The expansion of the Maxwell field can be p erformed in 
the same way as the metric. Although our results seem different with the EFT calculation in Ref. [16] . this is just 
the difference in the parameterization. We also confirmed the first law and the Smarr formula, which were shown for 
a sequence of charged black objects in Kaluza-Klein spacetime in Ref. [20| . 

The another way to construct the charged solution is using of the Harrison transformation pol l2l| which produces 
the charged dilatonic solution from the neutral seed solution. However, the charged stationary solution such like 
the rotating black ring cannot be produced by the Harrison transformation. In such case, one should rely on some 
perturbative methods with the Maxwell field as we used in this paper. Since the charged black ring is studied only 
by the blackfold approach [22[ , it is also interesting to study the charged black ring beyond the blackfold approach. 

We have not included the finite size effect which comes from the deformation of the horizon. To sec this effect, we 
must consider the multipolc perturbation which will be studied in the future work. 
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Appendix A: Gegenbauer polynomials 

The Gegenbauer polynomials C^x) are defined as the coefficients of the following generating function, 

— = y c?(x)t l , (ai) 
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where Cf(x) is written as 



It follows the equation 



(1 - x 2 )y" - (2i/ + l)y' + 1(1 + 2zy)7y = 0. 



(A2) 



(A3) 



Appendix B: Asymptotic limit of the Newtonian potential 

In the limit, r z,L, the Newtonian potential $ in Eq. (|21a[) behaves as 



<I> 



E 



A) 



2p% ^ L/r 



1 - 



2mLz 



Wo \ ^ 

r n-l (! 



Lr™" 1 ^ (1 + (mL/r) 2 ) 
2po f°° dt 



r 2 + m 2 L 2 r 2 + m 2 L 2 
s + 0(l/r n ) 



Lr n ~ 1 J (1 + i 2 )- 



^+G(l/r n ) 



(n — l)w n Lr T: 



r + 0(l/r"). 



1 



2mLz z 1 
r 2 + m 2 L 2 r 2 + m 2 L 2 



(r 2 + z 2 ) 1 - 



(Bl) 



Appendix C: Multipole perturbation in the near zone 

In this section, we give the equation for the multipole perturbation around the charged black hole. 



1. Ansatz 



The static spacetime can be written in this form, 

ds 2 = gtt{x)dt 2 + gij(x)dx l dx 3 . 



In this ansatz, we consider the linear perturbations around the black hole as = gp. v + h^, where the background 
metric is given by 



(CI) 



gf^dx^dx" = -f(p)dt 2 



dp 2 



P 2 dn 2 n+1 . 



(C2) 



Since the background metric has SO(n + 2) symmetry, we can expand the perturbations in terms of the spherical 
harmonics on the (n + l)-dimcnsional sphere. Here, we consider only scalar perturbations. We explicitly write all 
nonzero components of the perturbations, 



h pp — 



h a b 



1=0 

oo 

1=0 

oo 

l=i 

oo 1 oo 

p 2 V D,(p)[X> 2?6 - — - lab V 2 )Yt + p 2 lab V E^Yt, 
— n + 1 f— ' 



(C3a) 
(C3b) 
(C3c) 
(C3d) 



i=i 



1=0 
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where j a b is the metric of the (n + l)-dimensional unit sphere and D a denotes the covariant derivative with respect 
to 7 a b. Here, Y; is the spherical harmonics with Z-th multipole moment, which satisfies T> 2 Y[ = —1(1 + n)Yj. The 

gauge transformation for the metric perturbations is given by h^ v — > + C^g^v under the infinitesimal coordinate 
transformations — > x^ — Then, the components transform as 



Ai + , Bt -> B, + 2tf ' - Lfo d^d+ + fp%, 

p n + 1 



(C4) 



where we have expanded £ M with the harmonics Y/ and (i is defined as £ a = QT) a Yi. Note that we have not transformed 
the time coordinate t. Using the above gauge degrees of freedom, £ ; p and Q, we can set C; = D{ = Q. In the 
monopole case, Co and Dq are automatically zero and then £q becomes the residual gauge. For the vector potential, 
we write the perturbations, SA t , SA p and 5A a as 

oo 

6 A t = ^ a/ (p)Y ; , (C5a) 

oo 

SAp = J2Wp)Yh (C5b) 

1=0 

oo 

5A a = 5^c,(p)X> Y,. (C5c) 
i=i 

The U(l) gauge transformation of 5A^ — > SA^ + d^ip, where ip is a scalar function, transforms the potential as 

h -> b, + ip' t , ci^ci+ ipi (C6) 

where 4>i{p) is the expansion coefficient of ip for the spherical harmonics Yj. Therefore, we can set c; = (/ 7^ 0) 
by choosing ipi appropriately. In the monopole case, we also set 6 = 0, because Cq is automatically zero. The field 
strength becomes 



Z=0 

00 

*^ia = -£oiI> Y, I (C7b) 



00 



5F pa = -J2biV a Y l . (C7c) 



<5i? M; , = 8-7T 



<5TU - 

n 



2. Perturbation equation 

The Einstein equation for the linear perturbation is 

= SirSS^ (C8) 
Under the metric ansatz (|C3[) . the non- vanishing components of SR pL/ become 

SR tt = --L[-2p 2 fA'{ -{Zp 2 f + 2(n + l)pf)A' l +p 2 f'B' l 
-(n + l)p 2 /'^' - 2(p 2 /" + (n + l)p/')(4 - B«) 

+2Z(Z + n)A ; ]Y ; (C9a) 
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SR PP = ^j[-2p 2 fA'{-2(n + l)p 2 fE' l '-3p 2 f'A' l 

+( P 2 f' + 2(n + l) P f)B[ - (n + l)(p 2 /' + 4p/)E{ 

+2Z(Z + n)£ z ]Y, (C9b) 



^ = 77*l- 2 Pf A i-2npfE' l -(pf'-2f)A 



1 

(p/' + 2nf)B l ]V a Y l (C9c) 



«U = ^[V/^f - (P 2 /' + 2(n + - pf(A[ - B[) 

+2(pf + n/)(Bj - Sj) + Z(Z + n)£,]Y z7ab 

-i[^i + B l + (n - 1)£?!]X> I>6Y, (C9d) 
where we write only the Z-th pole components. On the other hand, the non-vanishing components of SS^ becomes 

SS tt = -f^—(2£a' l +£ 2 B l )Y l , (ClOa) 

71+1 

1 n 

SS PP = -——(2£a' l +£ 2 A l )Y l , (ClOb) 

J 71 T" -L 

<$Sa6 = P 2 7ab— ^-t [-2£a{ + £ 2 (£ ( - 4 - B,)]Yi, (ClOc) 
n + 1 

8S pa = -£a x V a Y h (ClOd) 
SS ta = f£b l V a Y h (ClOe) 

where £ is defined by Eq. (flgjl . From <5i? ta = and E q. (|C10ej> . we see 6; = for Z 7^ 0. Moreover, a; (Z 7^ 0) 
is determined by using the Einstein equation from Eqs. (|C9cp and (|C10d|) . Therefore, we do not need to solve the 
Maxwell equation directly to determine a; (Z 7^ 0). 

3. Master Equations 

From the traceless part of Eqs. (|C9d[) and (jClOcj) . we obtain the following algebraic relation 

i4,+B, + (n-l)f7, = 0, (Cll) 

because [(n + \) r D a T>} } — 7 a b£> 2 ]Y; 7^ when Z > 1. Note that for Z = 0, 1 we do not have any constraint from the 
perturbation equations but we can require this equation by using the residual gauge (See Eq. I|38p). 
Now, we introduce the following master variables, 

Ai = # - nil, Bi = IT - Ei =11, (C12) 

which satisfy Eq. (jCllj) . Using the master variables, we obtain the following master equations from the perturbation 
equations 

(n + l)p 2 fW + (n + l)p/n' - (n + 1)1(1 + n)U = -21(1 + n)*, (C13a) 

p 2 ft>" + (2p 2 f + (3n + 17)pf)W + (2n 2 - 1(1 + »))* = 0. (C13b) 
Using the new dimensionless variable x = (p/po) n , Eqs. (|C13j) become 

fx 2 ^ + 2fx^- Pl ( Pl + l)n = 2 — Vl (vi + 1)* (C14a) 

ax z ax n + 1 
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and 

^S + ^£ +4 ^f + (2-*( P « + l))* = 

(C14b) 

where = i/n. Using f(x) = (x — x+)(x — X-)/x 2 , Eq. (|C14b[) becomes 

d 2 ^ 

(x - x+)(x - x-)— + (4b - 2.T+ - 2a;_)— + (2 - + 1))* = (C15) 
ax z ax 

This equation has the three singular points : x = .x_,a;+ and oo. Therefore, the equation for \I/ becomes the 
hypcrgcomctric differential equation. Using £ — (x — X—)/(x+ — £_), it becomes standard form, 



d 2 V . 

(C16) 



- ^)— + ( 2 - 4 0— + (Pi + 2)(pj - 1)* = 



This is the same with one for the neutral case @. Similarly, defining n = xli./{x+ — X-), Eq. (|C14ap becomes the 
hypcrgcomctric differential equation with a source term 



f(i-f)^+«( P « + i)n = ^±^tt + e)* I 

n + 1 

(C17) 



where we defined e = — x_). So, our multipole solution, in general, can be written by four hypergeometric 

functions and the characteristic solution for Eq. (|C17|) . In the neutral case, we do not need to solve the equation for 
II, because II is obtained directly from ^ using Eqs. (|C9c[) and 5R pa = 0. 



Appendix D: Gauge transformation into the harmonic coordinates 



To match the near solution with the asymptotic solution in the overlap region, we need the near solution written 
in the harmonic gauge. If we write the harmonic coordinates as x^ix^) = (t, Ph{p)i Xi ■ • •)> the harmonic condition 
(V 2 £ M = 0) becomes 



V 2 (p h cosx) =0. 



(Dl) 



In the above, V 2 is the Laplacian for the full near metric g pv = gjft + hp V . Since the monopole perturbation does 
not depend on x, Eq. (|D1|) up to 0{1/L n ) becomes 



1 - W n) ) d P \( 1 + l hin) ) - B )d p 



(n + l)(l + E )p n - 1 p h , 



where =A + B + (n+ 1)E is the trace of h$ . 

We expand ph as pu = p^ + p\ , where p^ is of 0(1/L n ). Then, the O(l/L ) equation becomes 

d p {p n+1 fd pP f) = {n+l)p^pf. 
Assuming ph — > p as p — > oo, the solution is given by 



Ph 



2n p r ' 



The 0{l/L n ) equation is 



d p (p n+1 fd pPh n) ) - (n + l) P n - VL n) = d p (p n+1 fB d pPh 0) ) 



lp n+1 d p h^fd p p[ 0) (n + Vp^EvpI 



(0) 



(D2) 



(D3) 



(D4) 



(D5) 
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Note that the 0(1/L n ) equation is trivial in the gauge choice (|49l) . Here we consider cases for general d,di a o- 
Assuming p^ — > as p — > oo , the solution becomes 



As a result, we have 



p[ n) = 9 ^ £l ^k + OiptVp 2 - 1 )- (D6) 



ph(p) = p - a^r^ " n ° 2 c ^ + oiplVp 2 -- 1 ). (D7) 



This equation gives the coordinate transformation from (t, p, x, ■ ■ ■ ) into (t, Ph{p), X> • • • )• I n the harmonic coordinates, 
the near zone metric and the gauge field become 

5 t ( r r) = + MP)) = - (l - 2d - (1 - 2d - d) ^|) + 0(l/i 2 ", l/p 2 ,"), (D8a) 

ff&S? = f|^) 2 r 1 (p)(i + s„(p)) = i + ^ + 4rfi + — -d)+diA 2 'M/pn (D8b) 

w^h ydphj n npl \ n J 

= P 2 (! + *>00) = P 2 , ( 1 + ^ + 4l f 1 + — - d) + 0(l/L^, I/pi-)) , (D8c) 

47 r) = — %r+r(i-«o-2d + 0(i/i 2 'MMO)- (D8d) 
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